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Abstract 

The response of a cross-correlation measurement to an isotropic stochastic gravitational-wave 
background depends on the observing geometry via the overlap reduction function. If one of the 
detectors being correlated is a resonant bar whose orientation can be changed, the response to 
stochastic gravitational waves can be modulated. I derive the general form of this modulation as a 
function of azimuth, both in the zero-frequency limit and at arbitrary frequencies. Comparisons are 
made between pairs of nearby detectors, such as LIGO Livingston- ALLEGRO, Virgo- AURIGA, 
Virgo-NAUTILUS, and EXPLORER-AURIGA, with which stochastic cross-correlation measure- 
ments are currently being performed, planned, or considered. 



*Electronic address: jtwhelan@loyno.edu 



1 



I. BACKGROUND AND MOTIVATION 



A. Stochastic Backgrounds and Cross-Correlation Measurements 

One of the gravitational wave (GW) sources targeted by the current generation of 
ground-based interferometric and resonant detectors is a stochastic GW background 
(SGWB), produced by an unresolved superposition of signals of astrophysical or cosmo- 
logical origin. 

HQS 

Direct measurements can be made of, or observational limits set on, a 
SGWB by looking for correlations in the outputs of a pair of GW detectors. 

The response of a pair of detectors to a SGWB depends in part on the observing geometry, 
in the form of the relative orientation and separation of the two detectors. If one or both of 
the detectors is a resonant bar detector, its orientation can conceivably be changed, thereby 
altering the observing geometry and the response of the cross-correlation experiment to a 
stochastic GW background. For example, the experimental setup of the ALLEGRO resonant 
bar detector j4| allows for rotation of the detector. This produces a modulation of in the 
stochastic GW response which is used in the cross-correlation measurements being conducted 
between ALLEGRO and the interferometer (IFO) at the LIGO Livingston Observatory 
(LLOj.flflflS 

In the case of LLO and ALLEGRO, where the separation between the detectors (about 
40 km) is small compared to the wavelength of gravitational waves to which both detectors 
are sensitive (about 300 km), the azimuth dependence of the response is very nearly the 
same as for colocated detectors: sinusoidal with a period of 180°. j5|. This work considers 
analytically the general dependence of the stochastic gravitational- wave response on the 
orientation of a bar detector, and applies the results to existing pairs of detectors. 

B. Mathematical Details 

If a SGWB is isotropic, unpolarized, Gaussian, and stationary, it is completely described 
by its power spectrum. It is conventional to express this spectrum in terms of the GW 
contribution to the cosmological parameter Q = p/p CT it'- 

n ff\ 1 dp GW f dp GW 

"gwI/J -r; — 7 "77 • 

p crit d In / p crit df 
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If ^1,2 (/) are the Fourier transforms of the detector outputs, the correlation in the presence 
of a SGWB of spectrum Q(f) should be 
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(HfTHf)) = f^*(/ - f'M\m G w(\f\) (1.2) 



where 7(/) is a geometrical factor 111], LLlj involving the location and orientation of the two 
detectors, which we describe in detail in Sectional This factor also appears in the sensitivity 
of a standard cross-correlation technique^, for example in the case of a background whose 
^Gw(f) is constant over the frequency band of interest: 

II. DEFINITION OF THE OVERLAP REDUCTION FUNCTION 



^sens ~ ( T I df— ^TTT ) 



The factor j(f) appearing in f)1.2|) and ()1.3|) is typically (see e.g., |2j) defined as 

7(/) = ^ E ffd 2 Q A F 1A (n)F 2A (n)e l2 ^^-^ (2.1) 

A=+,x J J 

where f\ and f 2 are the positions of the two detectors, and Fia(ji) and F 2A (n) are their 
beam patterns. We follow the lead of by factoring the response tensors (cf. Appendix IXjl 
out of the integral and writing 

l(f)=d lab T« b d (a,s)d 2 cd (2.2) 

where 

s = ;}"?, (2.3) 
In - r 2 | 

is the unit vector pointing from one detector to the other, and 

27rf In — fo\ , 
a = - 1 1 - (2.4) 

c 

By using the definition of the transverse, traceless projector in (|A3J) . we note that we can 
write 

rS(a, *) = ^ // d ^ pTTfl c b d (2-5) 
Note that this gives an alternate definition of the overlap reduction function which is man- 
ifestly independent of any polarization basis: 



l(f) 



dlab d 2 cd J- J J d 2 ttn P TTha c b d e ^/n-(f 2 -rl)/c (3 g) 



Written in the form (|2.5|) . it is clear that T^(a, s) is symmetric and traceless on both pairs 
of indices ({ab} and {cd}). With s as the only preferred direction, there are only three 
independent tensors which can be created with these properties: 

Tit = (2-7a) 

r 2 S(s) = P T : b f s f s g p T Z (2.7b) 

T&{$) = P T t f S e S f s g s h P T9 c h d (2.7c) 

(Previous derivations included two additional terms which were not traceless, the 

coefficients of which are of course zero.) The most general possible form is thus 
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rS(a,S) = X;^(«)r»S (2-8) 

n=l 

and hence 

7 (/) = Pl ( a )^ o6 4 a6 + p 2 {a)dl ab s b s c d T 2ac + p 3 («)^ a ^ c S d 4 cd (2.9) 



The complete functional forms of the coefficients are derived in [1CJ and corrected in 
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121 ] . and a slightly simplified derivation appears in Appendix [B] of the present paper. 
First, however, it is elucidating to note the behavior at a = 0, which corresponds either to 
the low-frequency limit or to the case where the two detectors are co-located. In that case, 
()2.5|) becomes 

r2(o, = f[d 2 n h p™% (2.io) 



47T _ 

which is manifestly independent of the separation direction s. That tells us that it can 
contain only the s-independent term T^. I.e., P2(0) = = P3(0), leaving 

r^(o,s)= Pl (o)p T ^ = |- JJd 2 n n p™% (2.ii) 

It's easy to solve for pi(0) by taking the trace of 1)2.11)) and using ()A5|) and ()A6|) to say 

r^(o, §) = pi(o) • 5 = A ff d 2 n h ■ 2 = 10 (2.12) 
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This means that p(0) = 2 and 

7 (0) = 2^4 afe (2.13) 

Note that this makes, for instance, the demonstration that the overlap reduction function 
for coincident, coaligned IFOs with perpendicular arms is unity, extremely simple. [Just 
substitute (|A8|) in for both response tensors.] 



III. DEPENDENCE OF THE OVERLAP REDUCTION FUNCTION ON THE 
ORIENTATION OF A BAR DETECTOR 



In this section we consider in detail the form of the overlap reduction function between 
a bar detector with its long axis along the unit vector u and 

1. an IFO with axes along the unit vectors x and y 

2. a bar detector with its long axis along the unit vector x 

3. a general detector with response tensor d a t, 



A. Zero-frequency limit 

The relatively simple form (|2.13|) of 7(0) makes it worthwhile to consider briefly the 
overlap reduction function in this limit, which can in general be written as 

27r/|rW2 U o (3.1) 

c 

and thus also applies to detectors located at the same site. 



1. Correlations with an interferometer 

For correlations between an IFO and a bar, we use (jA8|) and (jA9|) and note that in this 
case df 3 is already traceless, so 

7 (0) = 2dfd 2ab = -[(£ • uf - (y ■ uf] (3.2) 

For the case where the detectors lie in the same plane, the IFO's arms are perpendicular, 
and the bar makes an angle of 9 with the IFO's "x arm" and \ — 9 with its "y arm" , so that 



x ■ u = cos# (3.3a) 

y ■ u = sm9 (3.3b) 

we have the familiar result that 

7(0) = cos 29 (3.4) 



2. Correlations with another bar 



In the case of correlations between two bar detectors, the response tensors are 

df = x a x b (3.5a) 
df = u a u b (3.5b) 

so that from f|7TTTT|) 

7 (0) = 2 (dfd 2ab - ^dA = 2 f{x • uf - = cos 26 + l - (3.6) 

where 6 is the angle between the two bars, so that once again x ■ u = cos 6. 

Note that the maximum value of (|3.fi|) . which occurs when 9 = 0, opposed to the 

maximum of unity for correlations between an IFO and a bar (13.41) or two IFOs. For this 

fin 

reason, the overlap reduction function for two bars is sometimes (see, e.g., |3j,lll|) normalized 
with an additional factor of | relative to (j2.1J) so that the maximum at zero frequency is zero 
for any type of detector pair. The normalization used in this paper is preferred, however, 
for several reasons: 

• Including a detector-dependent normalization factor in the definition of 7(/) would 
mean formulas like (jl.3|) would have different forms for different types of detector pairs. 

• While the maximum value of 7(0) for two bars is | rather than 1, its minimum value 
is — I rather than —1 (which is the minimum value for either two IFOs or an IFO and 
a bar). So in all three cases the amplitude of the orientation-induced modulation is 1. 

The latter point can be physically understood in terms of a bar detector being "more omni- 
directional" than an IFO. For an IFO, there is one optimal propagation direction for which 
the GW response is a maximum, namely a wave propagating perpendicular to the IFO plane. 
Two IFOs with their arms parallel will both respond ideally to waves propagating in this 
direction. The same is true if we rotate one of the IFOs 90°, except that the complete cor- 
relation of the signals has become an anticorrelation, so 7 = — 1. On the other hand, a bar 
detector has a one-parameter family of "optimal" propagation directions, all perpendicular 
to the bar. Two parallel bars share the same family of optimal directions. However, if we 
rotate one bar 90°, the planes perpendicular to the bars no longer coincide, and there is only 
a single optimal propagation direction perpendicular to the bars for which they both respond 



optimally. So the anticorrelation in the perpendicular arrangement is not as efficient as the 
correlation in the parallel arrangement, and |7(0)| is smaller for perpendicular bars than for 
parallel ones. 

B. General form 

Returning to consideration of the overlap reduction function 7(/) for general frequencies 
and detector separations, what interests us here is the dependence of j(f) on the orientation 
of the bar. We assume that the locations of the two detectors and the geometry of the other 
detector are fixed. The one variable is the azimuth of the bar, which is conventionally defined 
as an angle measured clockwise from the local geographic North. If we define unit vectors N 
and E pointing North and East, respectively, the orientation vector for a bar with azimuth 
C is 

u = iV cos C + ^ sin C (3.7) 

its response tensor is 

df = (N a cos C + E a sin C) (N b cos ( + E b sin C) = df + d% cos 2( + df sin 2( (3.8) 
where 

< 

4 
4 

This decomposition of the response tensor into a piece independent of the azimuth, plus 
pieces proportional to the sine and cosine of twice the azimuth, allows us to write the 
overlap reduction function as 

7(/) = 7o(/) + lc{f) cos2C + 7 s(/) sin2C (3.10) 

where 

lo,c,s(f) = d lab T^(a,s)d c d c ^ s (3.11) 

Given the location of a bar, it is relatively straightforward to work out the components in 
some convenient Cartesian basis of the unit vectors N and E at the location of the bar, and 
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N a N b + E a E b 
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N a N b - 


E a E b 


2 




N a E b + 


E a N b 


2 



(3.9a) 
(3.9b) 
(3.9c) 



thus construct the tensors , d^, and dg. Those can be used, along with the information 
about the geometry of the other detector and the separation between the two detectors, to 
construct the coefficients 7o(/), 7c(/), and js(f)- 

One final improvement on the form (J3.10j) is to define 7a(/) an d Cmax(/) according to 

7c(/) = 7A(/)cos2C max (/) (3.12a) 

75(/) = 7A(/)sin2C ma x(/) (3.12b) 

so that 

7(/) =7o(/)+7A(/)cos2(C-Cmax(/)) (3.13) 

At a given frequency /, 7o(/) is the orientation-independent piece of the overlap reduction 
function, 7a(/) is the amplitude of the orient at ion- dependent modulation, and Cmax(/) is 
the orientation (modulo 7r) for which the overlap reduction function is a maximum. Note 
that this is not necessarily the optimal alignment for stochastic background observations; if 
7o(/) < at the frequency of interest, the stochastic sensitivity is maximized by setting the 
azimuth to Cmax(/) + f (again modulo n). 



C. Examples for real- world detectors 



I3.13jl to pairs of real-world detectors. In 
1311 was used to calculate the orientation- 



This section contains several applications of 
each case the detgeom suite of matlab routines 
independent offset 7o(/), the amplitude 7a(/) of the modulation and the azimuth ( max (f) 
of maximum overlap, SIS db function of frequency. 

The major bar detectors around the world were oriented roughly parallel to one another 
as part of agreements through the International Gravitational Event Collaboration (IGEC) . 
This orientation is referred to here as "the IGEC orientation" . 

Fig. ^ shows the modulation for the ALLEGRO bar detector (Baton Rouge, LA, USA) 
correlated with the LIGO Livingston interferometer (Livingston, LA, USA). This is the 
closest pair of detector sites, separated by only 40 km. 

Fig. El examines the modulation for the AURIGA bar (Legnaro, Italy) correlated with 
the Virgo interferometer (Cascina, Italy). AURIGA is the closest bar detector to Virgo, 
separated by 220 km. 
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Overlap Modulation for LLO/ALLEGRO 
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FIG. 1: Modulation of the overlap reduction function for the ALLEGRO bar correlated with the 
LIGO Livingston interferometer. In the plot at left, the dashed line shows the amplitude 7a(/) 
of the azimuth-dependent oscillation, while the solid line indicates the azimuth-independent com- 
ponent 7o(/)- The detector sites are only separated by 40 km, so the amplitude and offset of the 
modulation change little between Hz and 1000 Hz. The sensitive frequencies for the correlation 
measurement are near 900 Hz, and the range of overlap reduction function values as the orientation 
is changed is —0.90 < 7(900 Hz) < 0.96, which makes this pair of detectors well suited to the mod- 
ulation of stochastic GW response. Likewise, the optimal orientation differs from that determined 
at zero frequency by only a fraction of a degree. The IGEC orientation, not shown on these axes, 
is -40, i.e., 40° West of North. 

Fig- El shows the modulation for the NAUTILUS bar (Frascati, Italy) correlated with the 
Virgo interferometer. The two sites are separated by 270 km. 

The result f)3.L3jl can also be applied to a pair of bar detectors. Fig. 0] shows the modula- 
tion for the AURIGA bar detector correlated with the EXPLORER bar (Geneva, Switzer- 
land). The detectors are separated by 460 km. EXPLORER is assumed to be in its actual 
IGEC orientation of 39° East of North, and the modulation is that associated with the 
AURIGA azimuth. 
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Overlap Modulation for Virgo/AURIGA 
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FIG. 2: Modulation of the overlap reduction function for the AURIGA bar correlated with the Virgo 
interferometer. The dashed line on the right-hand plot indicates AURIGA's IGEC orientation, for 
reference. The detector sites are separated by 220 km, so the frequency dependence of the overlap 
reduction function is important. It happens that in the azimuthal modulation 7a(/) is rather small 
at the frequencies of the bar's sensitivity (around 900 Hz). However, the azimuth-independent offset 
is non-negligible, which means that 0.14 < 7(900 Hz) < 0.22 for all orientations. In particular, it 
is not possible to change the sign of the gravitational-wave response by altering the orientation of 
a detector at the AURIGA site. 

IV. CONCLUSIONS 

For stochastic background measurements correlating a resonant bar detector with an- 
other GW detector, we have derived several analytical results about the dependence of the 
overlap reduction function j(f), and hence the response to an isotropic SGWB, on the bar's 
orientation. In particular: 

At zero frequency, the overlap reduction function for any pair of detectors is twice the 
contraction of the traceless parts of the detectors' response tensors (|2.13|) . For correlations 
between an interferometer and a bar lying in the same plane, this has a sinusoidal dependence 
on the angle between the bar and the interferometer's "x arm" (|3.4|) . For correlations 
between two bars, there is still a sinusoidal dependence on the angle between the two bars 
(13. 6|) . but now that modulation is offset by |, so that cross-correlation measurements with 
parallel bars 7(0) = | are more efficient than those with perpendicular bars 7(0) = — |. 
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Overlap Modulation for Virgo/NAUTILUS 



Overlap Modulation for Virgo/NAUTILUS 
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FIG. 3: Modulation of the overlap reduction function for the NAUTILUS bar correlated with 
the Virgo interferometer. The dashed line on the right-hand plot indicates NAUTILUS'S IGEC 
orientation, for reference. Since the detectors are separated by 270 km, the amplitude 7a(/) of 
the modulation has its first minimum below 900 Hz, and the range of overlap reduction function 
values is —0.18 < 7(900 Hz) < 0.07, depending on the bar's azimuth angle. Note also that while 
Cmax(900Hz) = 98°, this is the azimuth for which 7(900 Hz) is a maximum; the maximum of 
|7(900 Hz) | is in the perpendicular alignment, an azimuth of 8° East of North. Note also that 
while the actual IGEC azimuth of 44° is 36° away from this optimal orientation, this is not as 
close to a "null" alignment as it seems because of the offset 7^(900 Hz) = —0.06, which makes 
7(900 Hz) = -0.10 for the IGEC orientation. 

This can be understood, as explained in Sec. 1111 A 21 in terms of the distribution of optimal 
propagation directions for bar detectors. 

At higher frequencies, we have confirmed the empirical result seen in of a sinusoidal 
azimuthal modulation of the overlap reduction function, generally offset from zero mean 
()3.13|) . We have also verified that the modulation has a period of 180° with no higher-order 
moments, and give explicit expressions for the amplitude 7a(/), offset, 7o(/) and phase 
(max(f) °f the azimuthal modulation, in terms of geometrical quantities such as the detector 
response tensor and the unit vectors associated with the bar detector's location. 

Finally, we have illustrated these results in Section (jHI Cj) by applying them to several 
pairs of real-world detectors. LLO-ALLEGRO has the best observing geometry and the 
greatest opportunity for modulation, with 7(900 Hz) varying between —0.90 and 0.96 with 
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FIG. 4: Modulation of the overlap reduction function for the AURIGA bar correlated with the 
EXPLORER bar. EXPLORER is assumed to be in its actual IGEC orientation of 39° East of 
North. Note that, as discussed in Sec. 1111 A 21 the offset of the modulation at zero frequency is 
1/3 rather than zero. Note also that at zero frequency, the azimuth Cmax(O) of maximum overlap 
is nearly equal to the actual IGEC azimuth, which is to be expected since in that orientation the 
bars are nearly parallel. 

azimuth. The AURIGA bar has a reasonable observing geometry with Virgo, but little 
opportunity for modulation if it were rotated, with 7(900 Hz) ranging from 0.14 to 0.22, 
while NAUTILUS and Virgo have a slightly less favorable geometry on average, but greater 
potential modulation, with 7(900 Hz) between —0.18 and 0.07. 
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APPENDIX A: CONVENTIONS AND DEFINITIONS 

Lowercase Latin letters (a, b,c, . . .) will represent spatial indices ranging from 1 to 3, 
which can be raised or lowered with the Kronecker delta 5 a b- The Einstein summation 
convention will be used for these indices. 

The propagation direction of a plane wave is defined by a unit vector n; given such a 
unit vector, one can choose an arbitrary unit vector £ perpendicular to n and then define 
rh = n x £ so that {I, rh, h} form a right-handed orthonormal basis. Given such a basis, it 
is conventional to define transverse, traceless polarization tensors 

ef(n) = tt - m a m b (Ala) 
e a x b (n) = tm b + m a t (Alb) 

Note that while these tensors span the space of traceless, symmetric tensors transverse to 
n, they are not orthonormal under the obvious choice of inner product on that space, since 

e A ab(n)e%(n)=2 5 AA , (A2) 

where A and A' are polarization indices, which can take on the values + and x . 

It is useful to define the projector onto traceless, symmetric tensors transverse to the unit 
vector n: 

P™? d = \ E ^)e Acd (n) (A3) 

A=+,x 

The factor of 1/2 is needed to make this a projector, because of the factor of two appearing 
in (jA2|) . Note that while the polarization basis tensors depend on the (arbitrary) choice of 
£, the projector only depends on the unit vector h. Note also that, because it is a projector 
onto a two-dimensional subspace, 

pTTnab = 3 (A4) 

Another useful projector is that onto traceless, symmetric tensors: 

= Sf c S b d ) - l* a %<i (A5) 
Note that because this is a projector onto a five- dimensional subspace, 

= 5 (A6) 
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Finally, we note here the forms of the response tensors for interferometric and resonant 
bar detectors. The response of a gravitational wave detector to a GW-induced metric per- 
turbation h a b(t,r) is 

h(t) = h ab (t,r dct )d ab (A7) 

where r^et is the spatial location of the detector and d ab is a symmetric tensor describing 
the detector geometry. The response tensor for an IFO with arms along the (not necessarily 
perpendicular) unit vectors x and y is 

< fo ) = \{^ b - y a y b ) (as) 

while that for a cylindrical resonant bar detector with its long axis along the unit vector u 
is 

dlL) = n a u b (A9) 

Note that we choose not to make the response tensor for the bar detector explicitly traceless. 
It is often useful to work with the traceless piece of a response tensor 

d Tab = pTab jcd = jab _ 1 fiabjc (AlO) 

which for a bar is 

dj h % = u a u b - \8 ab (All) 

(d? b {o s is, of course, already traceless) 

Finally, a standard pair of quantities used in describing GW detector response is a de- 
tector's beam pattern functions for the two polarization states, defined by 

F A (n) = d ab e Aab (n) (A12) 

APPENDIX B: DERIVATION OF THE COEFFICIENTS IN THE OVERLAP 
REDUCTION FUNCTION 

In this appendix we find the values of the coefficients pi(a), P2(«), and /^(a) appearing 

in 

I3(M) = ^JJ d^P™? de ^ = j2Pn(a)T n ? d (Bl) 
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[cf. (12. 5|) and (|2.8|) ] We do this by contracting each of the {T n ^} in turn with T 



cd i 



a, si: 



TiX(a, $) = ^Jj d 2 n A P TT Zl e M (B2a) 

T 2 X(a, *) = ^ // d 2 n h P TT Z%s c e M (B2b) 

T s SrS(a, S) = A || P™« a <M c S d e**' (B2c) 
The three contractions appearing in the integrands can be written 

pTTnafe = g (gg^ 

P™s b s c = 1 - (n • s) 2 (B3b) 

P™3 w 6 *" = 5 [(e a + b (n) *A) 2 + « 6 (n) s a 5 b ) 2 ] = ~ [l - (n • sf] 2 (B3c) 
Where we have used 

^™ = ~ [e?(n) e +ac (n) + <"(n) e x ac (n)] = £ 6 4 + m 6 m c = 5 C 6 - n a n c (B4) 
and, defining a spherical coordinate system such that 

s = sin 9 cos £ + sin 9 sin rh + cos # n , (B5) 



e+ b (n) s a s b = (£ ■ sf - (rh ■ sf = sin 2 #(cos 2 - sin 2 0) = sin 2 fl cos 20 (B6) 

and 

e^(fl) s a Sb = s ■ e X n ■ s = 2(1 ■ s)(rh ■ s) = 2 sin 2 6* cos sin = sin 2 #sin 20 (B7) 
This makes the contractions 

TiK(a, s) = 5 f dfi = 10— = 10j„(a) (B8a) 

MM) = | Ai- = 5 + H*£) = 10 Mf!) (B8 b) 

2 ,/_! \ ar or ) a 

aftr«i/„^_ 5 [\i -o,,!,,,^^.^^ 4sina 12cos« , 12 sina^ _ on j 2 (a) 



rs25rs(a, s) = - / (1 - 2/ + //) e ^ ^ = 5 -— 5 -j— + — g- = 20 , 

(B8c) 

To work out the contractions of the {T n ^} with the right-hand side of ()B1|) . we just need 
to contract each of them with the others. We demonstrate the calculation explicitly here in 
our notation: 

T^K = P T H = 5 (B9) 
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To get 



we note that 



so 



Next 



To work out 



We note that 



so 



Ti abrp cd ryVab £ ~c 
l C d 1 2 a b - r ac S bS 



_S a 5 b + -5 a 5 b ' 
2 a c 1 2 c a 3 



i 5 ™ = ^« + M - -%5 ab 5 ac =[^ + 1-1)6 



3 1 1 
2 + 2 ~ 3 



lcd 1 2ab 



rp abrp cd pTabg g gcgd i 

I lcd I 3ab - r cd S a SbS S - 1 - 



1 _ 2 

3 ~ 3 



rp abrp cd_ pTabg/g -eg pTep 
1 2 c d 1 2 a b — r ef S S b S ac 



rp abrp cd 
1 2 c d 1 2ab 



1 1 1 

' ! +3 + 3 + < 



To calculate 



we note that 



so 



Finally, 



rp abrp cd _ gagbpTe/g ggpTcdgcgd 
I 2 cd J-3ab - S S r ab S f S f eg S S 



Ti abrp cd 
2cd I 3ab 



Ti abrp cd 
3cd 1 3 a b - 



Sah - ^S ab ) I S a S b - -5 ab 



1 1 1 3 
X -3-3 + 9 



We can summarize these results as 



^ rp ab^ 

Ti ab 
2cd 

\Ts? d J 



rp ab rp ab rp ab 
1 lcd 1 2cd 1 3cd 



' 5 5/3 2/3^ 
5/3 17/18 4/9 
^2/3 4/9 4/9^ 
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which we can use, along with 



r~p ab T ab r~p ab 
I lcd I 2 c d I 3 c d 



to combine the three sets of contractions in the matrix equation 



^10j («^ 



10 
20 



h (g) 

a 



T ab^ 
1 lcd 

T< ab 
2cd 

rp ab 



( 



abK a i s ) 



5/3 2/3 



5/3 17/18 4/9 
^2/3 4/9 4/9 J 



Inverting the matrix gives 



/02(a) 



/ 1 



: 1 \ 

•1 4 -- 

5 35 

2 8 / 



20^ 



/ 



5 -10 5 
-10 40 50 



5 



which are the standard coefficients in the expansion ()2.9)1 



i 9 J 



25 



/02(a) 

^s(a) y 

a 

32(g) 
\ a* / 



(B22) 



(B23) 



(B24) 
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